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l/^ i Abstract It's generally believed that young and rapidly rotating pulsars are important sites 

QQ ■ of particle's acceleration, in which protons can be accelerated to relativistic energy above the 

" ■ I polar cap region if the magnetic moment is antiparallel to the spin axis(/x • Jl < 0). To obtain 

f^*) , the galactic diffusive neutrinos and gamma-rays for TeV, firstly, we use Monte Carlo(MC) 

method to generate a sample of young pulsars with ages less than 10° yrs in our galaxy ; 

secondly, the neutrinos and high-energy gamma-rays can be produced through photomeson 
^^ ' process with the interaction of energetic protons and soft X-ray photons (p + 7 ^- A+ ^■ 

j^ I n+7r+/p+7r") for single pulsar, and these X-ray photons come from the neutron star surface. 

The results suggest that the diffusive TeV flux of neutrinos are lower than background flux, 

which indicated it is difficult to be detected by the current neutrino telescopes. 
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1 INTRODUCTION 

Very high-energy [VHE,ii^ > lOOGeT^] neutrinos and gamma-rays from astrophysical objects can provide 
a clear indication of the origins of galactic and extragalactic cosmic rays. The radiations probably come 
from gamma-ray bursts, active galactic nuclei,pulsars etc. In our Galaxy most high-energy gamma-rays and 
neutrinos are associated with pulsars, supernova remnant and pulsar nebula(Bednarek et al. 2005;KistIer 
& Beacom 2006). Young pulsars are persistent and periodic sources; they also have shorter distances to 
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searching for the signals. If the TeV neutrinos and gamma-rays can be detected from pulsars, it will help 
improve our knowledge about the hadronic process taking place in the magnetosphere. 

Recently, some studies (Link & Burgio 2005, 2006;Bhadra & Dey 2009;Jiang et al. 2007) have shown 
that if pulsar(age less than lO^yr) magnetic moment is antiparallel to the spin axis (fift < 0),A-resonance 
will be produced when pulsars-accelerated ions interact with these thermal radiation field. Then high-energy 
neutrinos and gamma-rays decay by A-resonance. Link&Burgio(hereafter LB) considered that near the 
surface of the neutron stars, the protons or heavier ions can be accelerated by polar caps to PeV energy. When 
the thermal radiation flied of pulsar interact with accelerated ions with PeV energy, A-resonance state 
may occur. This process is effective, and pions subsequently decay to muon neutrinos or gamma-rays . LB 
calculated energy spectrum of muon neutrinos from some pulsars and estimated the event rates at Earth. The 
spectrum sharp rise at ^50 TeV, corresponding to the onset of the resonance. The flux drops with neutrino 
energy as e^^up to an upper energy cut-off that is determined by either kinematics or the maximum energy 
to which protons are accelerated. Bhadra & Dey(2009)estimated TeV gamma-ray flux at the Earth from a 
few nearby young pulsars and compared with the observation, they found that proper consideration of the 
effect of the polar cap geometry in flux calculation is important, then they revised the event rates. 

In this paper we find the rj (ratio of polar cap area to neutron star surface area)is double the Bhadra & Dey 
give. We use revised polar cap geometry to calculate TeV neutrinos and gamma-rays spectrum from some 
young pulsars, and estimate the event rates at Earth. In order to estimate diffuse TeV neutrinos and gamma- 
rays radiation from Galactic Plane, a galaxy pulsar sample is required. It is the better way to simulate 
a pulsar sample with MC method (Cheng & Zhang 1998;Zhang & Harding 2000;Jiang et al. 2007). At 
last,we obtain the TeV diffuse Neutrinos and Gamma-rays flux from young pulsars in our galaxy. 

In section 2, we reviewed the model of LB. In section 3,we calculate neutrino and gamma-ray spectrum 
from single pulsar by revised polar cap geometry.In section 4 and 5, we obtain a sample of young pulsars in 
our galaxy with Monte Carlos method, then estimate the diffuse neutrinos and gamma-rays emission from 
Galactic young pulsars. At last the discussion and conclusion are presented. 

2 THE MODEL OF NEUTRINOS AND GAMMA-RAYS FROM YOUNG PULSARS 

Neutron stars have enormous magnetic fields(B > lO^^G) and high rotation rates(tens of Hertz)which act as 
a very powerful generator. Charges will be stripped off the highly conductive surface and accelerated some- 
where above the stellar surface. The acceleration of particles mechanisms from pulsars generally has been 
divided into the polar gap model (Ruderman & Sutherland 1975;Arons & Scharlemann 1979;Daugherty 
& Harding 1996;Gonthier et al. 2002) and outer-gap model(Zhang & cheng 2004;Cheng et al. 1986). In 
the former, particles are accelerated in charge-depleted zone region near the magnetic pole of the neutron 
star. In the outer-gap, it will take place in the vacuum gaps between the neutral fine and the last open line 
in the magnetosphere. Therefore, acceleration region in the polar-gap model is close to the neutron star 
surface, whereas in the outer-gap model it is near to the light cylinder 

In polar gap model, because of existing large rotation-induced electric fields, particles can be ex- 
tracted from the polar cap surface, and then be accelerated, finally form the primary beam. The po- 
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A(/) = BgR'^fl^ /2c^{Goldieich & Julian 1969),where Bg is the strength of the dipole component of the 
field at the magnetic poles {Eg ^ lO^^G) and R = 10^i?6 is the radius of the neutron star and p,„ is the 
spin period in milliseconds ,51 = 2tt/p is the angular velocity(wherep is the period ) and c is the speed of the 
light. The magnitude of the potential drop A0 could be 7 x lO^^Bi2P,:n^ volts (Bs = B12 x 10^^)(Goldreich 
& Julian 1969). L B consider that If electric filed is no or little screening and /i ■ 17 < ( expected to hold 
for half of total pulsars), ions and protons will be accelerated to Pe V energy near the surface of pulsars, /x is 
magnetic moment, O is angular velocity. Protons or ions accelerated by pulsars will interact with the ther- 
mal radiation of pulsars. If they have sufficient energy and exceed the threshold energy for the A-resonance 
state(A+ is an excited state of proton, with a mass of 1232MeV).The A-resonance state may occurr. The 
threshold condition for the production of A-resonance state in p — 7 interaction is given by 

epe^(l - cosOp.,) > O.SGeV^, (1) 

where ep and e^ are the proton and photon energy, and 9p^ is the incidence angle between the proton 
and the photon in the lab frame. Young pulsars have typically temperatures of Too — O.lkeV, and photo 
energy e^ ~ 2.8fcToc(l + Zg) ^ OAkeV, where Zg w 0.4 is the gravitational red shift and Toois the 
surface temperature measured at infinity. In a young pulsar's atmosphere, the condition of production A- 
resonance is -Bi2^ms^o.ifeeV > 3 x 10^^(Link & Burgio 2005;Link & Burgio 2006), where To.ifeeV = 
{kToo/ O.lkeV), Too ^ O.lkeV is typical surface temperature of young pulsars. This condition holds for 
many young pulsars, so A-resonance could existed in many pulsar's atmosphere. Gamma-rays and neutrinos 
subsequently decayed through following channels 

p + 7 -^ A+ ^ { 

^7r+ — !• 71 + e+ + i^e + fp + i7^. 

In this part, based on the L B model;we estimated emitting neutrinos and gamma-rays from pulsars. The 

flux of protons accelerated by polar gap can be estimated 

Ipc = cfd{l - fd)noApc, (2) 

where no{r) = BsR^Q,/{ATiZecr^) is the Goldreich- Julian density of ions in distance r,fd < 1 is the 
fraction depletion in the space charge in the acceleration region. It's a model dependent quantity {fd = 
correspond to no depletion and /^ = 1 is full depletion), so the density in the depleted gap can be written as 
/(i(l — /d)7io,where Ape = ri2TiR^ is the polar cap area, rj is theratioof polar cap area to the half of neutron 
star surface area, when ?] = 1 is the half of neutron star surface area. In LB(2006)and Jiang(2007)they use 
hemisphere surface area to calculate. The typical radius rpc = R{riR/cY/'^ (beskin et al. 1993),the polar 
gap surface is Ape ~ t^i'^c — 7rili?'^/c,so ry = Q,R/{2c). So the maybe 77 = ili?/(2c) is more appropriate. 
For young pulsars with surface temperature T^o, the photon density close to the neutron star surface is 
n^{R) = (a/2.8fc)[(l + Zg)roo], ais the Stefan-Boltzman constant. Numerically 7i-y(i?) c^'^xlO^'^T^^^^^. 
At radial distance r, photon density will be n.y(r) = n.y{R){R/r)'^. The probability that a PeV energy 
proton staring from the pulsar surface will produce A particle by interacting thermal filed which is given 
by(link&Burgio 2005) Pc ~ I — JrP{''')- Where dP/P = —n^{r)ap^dr. Thus the total flux of gamma 
ray/neutri generated in pulsars from the A+ resonance is 
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Fig. 1 The neutrino flux displayed for crab,vela, the case of Unear(solid line ) and quadratic 
proton acceleration (doted line). 



where ^ is 4/3 and 2/3 for gamma rays and muon neutrinos,respectively. At the distance of d, the phase 
averaged gamma ray/neutrino flux at Earth from a pulsar is 



2c^Cv.fbfd{i - fdho ( I ) Pc, 



(4) 



where ft is duty cycle of neutrino or gamma-ray beam, <^ is the effect due to neutrino oscillation (the decays 
of pions and their muon daughters result in initial flavor ratios 0^_. : (f>i, : 0^^ of nearly 1:2:0 but at 
large distance from the source the flavor ratios is expected to become 1:1:1 due to maximal mixing of i^^ 
and !/,-). C = 1 and 1/2 is gamma rays and muon neutrinos respectively. We want to obtain the spectrum of 
neutrinos and gamma-rays, using the following differential from 



-— ^ 2c^Cvfbfd{l - fd)no -J -^— , 



(5) 



Taking fd = 1/2 and Z = A = 1 for estimating upper limits on the flux. The neutrino and Gamma-rays 
energy flux are estimated as 



dP 

^ = 3 X 10-»^Cr^f,Bi2P;nld]:lTo.ikev-T^. 



(6) 



de, .....^..^^."fcpc— - ^^ 

Where dkpc is the distance of pulsar to earth, and dPc/dx is the probabihty of proton converting to A+ per 
unit energy interval; the details have been shown in equation(26) at Link & Burgio(2006). 

3 TEV NEUTRINOS AND GAMMA-RAYS SPECTRUM 

We use the equation (6) to estimate flux of neutrinos and gamma-rays and take Z=l and f^ = 1/2 
throughout this work. Taking linear (7 = 1) and quadratic (7 = 2) proton acceleration laws. Linear 
acceleration is corresponding to an accelerating filed which is constant space and quadratic accelerating 
field grows linearly with height above the star. We calculate the neutrinos and gamma-rays energy flux of 
crab,vela,PSRB 1509-58 andPSRB1706-44,when7/ = Oi?/2c showed in the fig. 1 andfig.2. The parameter 
of source presents in Table. 1 . 



For either acceleration law, the spectrum turns sharply at e^, 
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Fig. 2 The gamma-ray flux displayed for crab,vela,the case of linear(solid line) and quadratic 
proton acceleration (doted line). 

Table 1 The parameter of some pulsars 



Source 



kpc 



P 
ms 



Bi2 

G 



Tb.lfce 



h ZATt (LB06) 



dAdt 

km~ yr 



dN 
dAdt 



km yr 



Crab 


2 


33 


3.8 


-1.7 


0.14 


45 


0.151 


Vela 


0.29 


89 


3.4 


0.6 


0.04 


25 


0.0313 


B1706 - 44 


1.8 


102 


3.1 


1 


0.13 


1 


0.0038 


B1509 - 58 


4.4 


151 


0.26 


1 


0.26 


5 


0.0057 



Table 2 The integral Tev gamma-ray fluxes comparison between 
the predicated and the observed. The observed upper limits for 
Crab, Vela,PSR B 1509-58 and PSR B 1706-44 (Aharonian et al. 
2006) 
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conversion becomes restricted. The flux is lower about 3 magnitude when we use pulsar surface to calculate 
.But it is more consistent with the observed upper limits of gamma-ray fluxes. numerical values of the 
integral TeV gamma-ray fluxes are obtained for pulsars listed in Table 2. 

In table 2, The gamma-rays flux estimated with ry = lis obviously higher than the observed upper limits. 
But with rj ~ ilR/2c it is a little lower and more consistent with the observed upper limits; Thus taking 
1] = VtR/{2c) is more reasonable. 

Large-area neutrino detectors use the Earth or ice as a medium for conversion of muon neutrino to 
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interactions below a detector on the surface of earth. We can use the flux of neutrino estimate the count rate 
in the detector The conversion probability in the Earth is(Gaisser et al. 1995). 

P,^^^ = 1.3 X lO-'\eJlTeV), (7) 

The muon event rate is 

-dXdt^ J '^'"ir/"'^^^- ^^^ 

In Table 1, expected count rates what we estimated present in the last column of for the crab,vela,PSR 
B1509-58 and PSR B1706-44,use L = 0.1 and linear acceleration. Expected count rates (LB06)are shown 
in the penultimate column. Compared with the expected count rates shown in the last column, expected 
count rates(LB06)are higher It means that the detectable by IceCube will not look bright. (Abbasi et al. 
2012) It shows that IceCube data severely constrain these optimistic predictions of LB, and (Bhadra & Dey 
2009) pulsars are unlikely to be strong source of TeV neutrino. So proper consideration of the effect of polar 
cap geometry in flux calculation is important. 

4 THE DIFFUSE NEUTRINOS AND GAMMA-RAYS EMISSION FROM GALACTIC YOUNG 
PULSARS 

Over 1 800 radio pulsars are known (Manchester et al 2005), more than 400 pulsars ages are less than lO^yr, 
most of which are candidates. If the magnetic moment is antiparallel to the rotating axis(yu • Vi < O),ions can 
be accelerated in the charge-depleted gap near the star surface. So young pulsars are also potential neutrino 
sources, although some of them are likely weak source, the total contributions maybe significant. Therefore 
we use MC method simulate pulsars sample with ages less than 10^?;r,then determine which is neutrino 
pulsars in sample. We also estimate the diffuse neutrino flux and gamma-ray flux. In other words,we will 
cumulate all potential neutrino pulsars in the sample. 

4.1 Pulsars Sample simulated by Monte Carlo method 

We produce the Galactic young pulsar population by following assumptions. (e.g.Sturner & Dermer 
1996;Cheng & Zhang 1998;Zhang & Harding 2000;Zhang & cheng 2004; Jiang & Zhang 2006). Pulsar 
sample is simulated by MC method which is described by following steps in this paper (see Jiang et al. 
2007 for detail). 

Following the steps, we use the conventional assumption for Galactic pulsars birth rate (^ 2/lOOyr) 
and enlarge 10 times, generate about 160000 pulsars in Galactic plane during past millions years. In this 
sample about 3890 pulsar could be detected radio. A lot of pulsars are not be observed because of selection 
effect. Compared with the ATNF catalogue pulsars, our simulated sample is consistent with the distribution 
of the observation. The normalized histogram distribution of pulsar period, pulsar surface magnetic field 
and distance are shown in fig. 3. The shadow histogram represents the ANTE sample,the solid histogram 
represents simulated sample. We also compare the distribution Galactic Longitude and Galactic Latitude 
with the observed sample which is shown in fig.4. The shadow histogram represents the ANTE sample,soUd 
histogram represents the simulated sample. Obviously,most pulsars distribute in Galactic Latitude \b\ < 5° 
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Fig. 3 Normalized histogram distribution of pulsar period, distance,surface magnetic field and 
period derivation. Shadow histogram represent observations sample sample. Solid histogram rep- 
resent simulated sample. 



4.2 The Diffuse Emission of Neutrinos and Gamma-ray from Young Pulsars 

The diffuse neutrinos and gamma-rays emission, both in Galactic and extragalactic are very interest for 
astrophysics, particle physics, and cosmology. The diffuse Galactic emission is produced by interac- 
tions of cosmic rays, mainly protons or electrons interact with the interstellar gas (via tt*' -production and 
bremsstrahlung) and radiation field (via inverse compton scattering). We estimate the gamma-rays and neu- 
trinos flux((i(/)„ i)/(ie!y by using equation(6). Then we can obtain the flux of neutrinos and gamma-rays from 
all pulsars 



N 






i=l 



dei,(Kli 



(9) 



where N is the number of gamma-ray or neutrinos pulsars, and dQ,i is the solid angle for the i pulsar, 
(Kli ~ '^irfb.i ^ 1- The results of neutrino energy flux are shown in fig. 5. The energy flux sharply increase 
at about 50Tey,corresponding to the onset of the resonance. Most energy is emitted between 50TeV and 
O.STeV. After the onset of the resonance,the spectrum drops approximately as e^^, because the phase space 
for conversion becomes restricted. Just hke the spectrum for a single pulsar,it is a typical of first-order Fermi 
acceleration(Abbasi et al. 2012). The linear proton acceleration is about 5 times than quadratic acceleration, 
so the linear proton acceleration is the main acceleration. From the figure 5, it will be not easier to search 



Z.-X. Li et al. 




-100 -50 o 50 100 

Galactic Longitude(Deg.) 




Galactic Latitude(Deg.) 



Fig. 4 Normalized histogram distribution of pulsars in Galactic Longitude and Galactic Latitude. 
Shadow histogram represents observations sample. Solid histogram represents simulated sample. 



little. For comparison, we also give the predicted fluxes with sensitivities of AMANDA-II, ANTARES, and 
IceCube detectors. AMANDA detector reached a sensitivity of 8.9 x lO^^GeVcm^^ s^^ sr^^ after the first 
4 years of operation in 2000-2003 (Hill 2006). The ANTARES detector sensitivity reached 7.8 ± 0.99 x 
lO^^GeVcm^^ s^^ sr^^ after 1 year of data on diffuse flux, and 3.9 ± 0.7 x 10~^GeVcm^^s^^sr^^ 
after 3 years of data on diffuse flux (Montaruli 2005). The IceCube detectors have reached a sensitivity of 
(2 — 7) X lO^^GeVcm~^ s~^ sr~^ after 3 years of operation (Ribordy 2006). It sees clearly that it not 
easier to detect the signals of neutrino flux from young pulsars, because the neutrino flux energy is lower 
than the background atmospheric neutrino flux. From Figure 5, we also compare the predicted diffuse 
neutrino fluxes from Active Galactic Nuclei (AGNs)(Stecker 2005)and Gamma Ray Bursts (GRBs)(Liu & 
Wang 2013). Whereas we can find the diffuse neutrino fluxes from AGNs are larger than young pulsars and 
GRBs at the energy from 60TeV to 3PeV. The diffuse neutrino fluxes from GRBs and young pulsars at the 
ranges from lOOTeV to IPeV almost the same and lower than ATM. The distribution of estimated diffuse 
neutrino flux from these young pulsars is narrower than AGNs and GRBs; the flux from AGNs and GRBs 
peaks at about 20PeV and IPeV respectively(He et al. 2012). The neutrino production in young pulsars is 
constrained by the accelerating efficiency of the protons, while the AGNs are more powerful accelerators 
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this energy range. The predicted diffuse neutrino from pulsars in Galaxy(Jiang et al. 2007) is easier to be 
detected, because the flux is above the sensitivity threshold of IceCube. Why they obtained higher flux than 
us,because the polar cap area we calculated is different. 



£ 
o 

> 

CD 



-e- 
■o 



O 




.-13 I ' i II 



4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0 6.2 6.4 6.6 6.1 

Log,„(B;GeV)) 



Fig. 5 The estimated diffuse neutrinos flux from young pulsars in our 
Galaxy. Solid and Dot line correspond to the fluxes computed for the 
cases of linear and quadratic proton acceleration. The solid lines labeled 
(l),(2),and(3)represent the sensitivities of AMANDA-II,ANTARES, and 
IceCube, respectively. The lines labeled (4),(5)represent the predicted fluxes 
from AGNs(Stecker 2005)and GRBs(Liu & Wang 2013) The region labeled 
"ATM" is the background atmospheric neutrino flux(Ribordy 2006). 

The results of diffuse gamma-rays flux are shown in fig. 6, Solid and Dashed line correspond to the fluxes 
computed for the cases of linear and quadratic proton acceleration. 



5 DISCUSSION AND CONCLUSION 

To summarize, based on the LB model and Bhadra(2009), we calculated the neutrinos and gamma-ray 
spectrum of some pulsars with appropriate polar cap area, we use MC method simulate pulsar sample with 
ages less than lO^yr and ji-H. < in our galaxy. Compared with the ATNF catalogue pulsars, our simulation 
is succeed. According to the sample, we also estimated the diffuse neutrinos and gamma-rays in our galaxy. 
Our results show that the diffuse flux distributing in the energy range from ^ 50TeV to ^ 3PeV is lower 
about 3 magnitude than Jiang(2007)only due to different methods on dealing polar cap area. 

These results show that the signals of neutrino flux from young pulsars are not easier to be detected, 
because the energy range about ^ bOTeV to ~ 3PeV is lower than the sensitivity threshold of IceCube 
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Fig. 6 The estimated diffuse gamma ray flux from young pulsars in our 
Galaxy. Solid and Dashed line correspond to the fluxes computed for the cases 
of Unear and quadratic proton acceleration. 

not been detected from the direction of any pulsar by the AMANDA-II telescope(Ackermann et al. 2005; 
Ackermann et al. 2008 ; Ahrens et al. 2004). So proper consideration of the effect of polar cap geometry in 
flux calculation is important and pulsars are unlikely to be strong source of TeV neutrinos. Compared with 
the TeV diffuse neutrinos flux emitted from AGNs(Stecker 2005) and GRBs(Liu & Wang 2013), AGNs 
maybe be the primary TeV neutrinos source. 

The energy data on the diffuse gamma-rays fluxes from Galactic plane around ~ bOTeV to '^ 3PeV 
have not been found now. So we cannot give the percentage of gamma-rays produced through A-resonance 
among the whole diffuse gamma-rays in our galaxy. We expect experiments like the proposed High Altitude 
Water Cherenkov (HAWC) detector will be constantly used to survey large regions of the sky, in particular 
the Galactic plane, at gamma-ray energies up to lOOTeV with 10 to 15 times the sensitivity of Milagro. If 
we have the energy data on the diffuse gamma-rays fluxes from the Galactic plane around ^ 50TeV to 
~ 3PeV. We can estimate the percentage of pulsar's contribution on this energy range. 

This work have some uncertainties, such as pulsar's birth rate and how many pulsars satisfy /i • O < 0. 
Estimating thermal photons from the neutron star surface is not well determined and the physical processes 
in the magnetosphere are still not clear, resulting in some uncertainties in the LB model as mentioned by 
(Link & Burgio 2005, 2006). Therefore, in order to improve the prediction of the diffuse muon neutrino 
and gamma-ray flux from young pulsars, it is crucial to understand more about the neutrino flux and spectra 
emitted by a Single pulsar. 
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